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I . Abstract 
5h . 

Starting from an action for discretized gravity we derive a canonical formalism that ex- 
actly reproduces the dynamics and (broken) symmetries of the covariant formalism. For 
£Nj i linearized Regge calculus on a flat background - which exhibits exact gauge symmetries - we 

derive local and first class constraints for arbitrary triangulated Cauchy surfaces. These 
constraints have a clear geometric interpretation and are a hrst step towards obtaining 
anomaly-free constraint algebras for canonical lattice gravity. Taking higher order dynamics 
into account the symmetries of the action are broken. This results in consistency conditions 
on the background gauge parameters arising from the lowest non-linear equations of motion. 
In the canonical framework the constraints to quadratic order turn out to depend on the 
| background gauge parameters and are therefore pseudo constraints. 

These considerations are important for connecting path integral and canonical quantiza- 
tions of gravity, in particular if one attempts a perturbative expansion. 

X 

H ! 1 Introduction 

Many approaches to quantum gravity, in particular path integral approaches such as Regge 
calculus [I], introduce an auxilary discretization as a regulator into the models. Although this 
facilitates the construction and finiteness of the models, these methods jeopardize diffeomor- 
phism symmetry, which is the continuum symmetry of general relativity. This symmetry is 
deeply entangled with the dynamics of the theory, therefore it might be fruitful to preserve a 
notion of diffeomorphism symmetry as far as possible. This could largely constrain possible 
quantizations of discrete models, moreover, it could provide a tool for controlling lattice effects, 
or, in other words, the independence from the chosen discretization (usually a triangulation), see 
for instance the discussion in [2]. Taking care of diffeomorphism symmetry can help to obtain 
the correct semi-classical limit and, in particular, to obtain the correct degrees of freedom in 
the large scale limit (as breaking of gauge symmetries introduces additional degrees of freedom). 

There has been some discussion in the literature [31 [U [5j [6] whether there exists an exact 
or only approximate notion of diffeomorphism invariance, or more generally gauge symmetry, 
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in discretized gravity, in particular in (the covariant form of) Regge calculus. For flat solutions 
it has been known [9] that symmetries exist: any vertex in the bulk of the triangulation 
can be translated in the embedding flat space without changing flatness, and therefore without 
leaving the space of solutions to the Regge equations. For a regular lattice these symmetries can 
be connected to the gauge modes in linearized gravity. 

But in [TO] it was shown with the help of an explicit example that these gauge symmetries are 
broken for Regge solutions with curvatureQ One might take the conclusion that generically in 
Regge calculus there are no symmetries and that the symmetries on flat space are an exception, 
not very relevant for the framework. As shown in [8j, however, the symmetries on the flat 
background are essential for obtaining the correct number of degrees of freedom in the continuum 
limit. Thus, assuming that generically there are no gauge symmetries and hence that all the 
degrees of freedom are physical would lead to a quite singular understanding of the continuum 
limit. Indeed, in toy examples one can show that the pseudo gauge modes behave dynamically 
very differently from the true physical modes, having a vanishing kinematical term Also 
in the Regge action expanded on a flat background the pseudo gauge modes obtain only a 
non-vanishing contribution from the (higher than second order) potential terms. 

We therefore rather speak of a symmetry breaking induced by discretization and, very im- 
portantly, by the choice of the discretized action. In fact, for 3d gravity with a cosmological 
constant the standard Regge action leads to symmetry breaking, whereas the so-called perfect 
action, corresponding to choosing simplices with homogeneous curvature instead of flat ones for 
the discretization, does not [TO] . 

The breaking of symmetries has particularly severe implications for the canonical formalism. 
Gauge symmetries of the action lead to constraints in the canonical formalism. These are the 
generators of the gauge transformations and hence form an algebra which is first class. The 
breaking of symmetries should therefore affect the constraints in some way. 

Indeed, performing the canonical analysis of the continuum action first and then discretizing 
the resulting constraints leads to additional terms that convert the continuum first class con- 
straint algebra into a second class constraint algebra! see for instance [HI E]. Note that this 
approach usually involves a change of set-up for Regge calculus, namely a continuous time, but 
still a discretized space [12] . An second class constraint algebra means that the constraints are 
not automatically preserved under time evolution. On the classical level one could deal with 
this issue by fixing the gauge parameters, that is lapse and shift, so that the constraints are 
preserved by time evolution |12|, I14j. The situation is, however, much more complicated in the 
quantum theory, see for instance (6j [15] and references therein. 

Nevertheless, there are attempts to derive a set of first class constraints for discretized 
theories. So far these succeeded only for 3d gravity, where the theory is topological, or for 
4d gravity in the sector where the triangulation is such that only flat space solutions arise 
[161 H7j . As the example of 3d gravity with a cosmological constant [TO] shows, the breaking 
of symmetries, and therefore the appearance of discretization anomalies, is not per se bound to 
discretization, but depends on the choice of the discrete action and, as a result, on the discrete 
dynamics. There might also exist a choice of discretized constraints in the canonical framework 
which are first class. In fact, in [TO] it was shown that once a discrete action with symmetries 
has been found one can derive first class constraints. 

To use this result, one has to develop a canonical formalism which exactly reproduces the 

More correctly, one expects the translation symmetry to be broken for vertices adjacent to triangles which 
have a non-vanishing deficit angle, that is curvature. But also in solutions with curvature there might exist 
vertices, for which all the deficit angles at the adjacent triangles vanish. For these the translation symmetries are 
preserved. 

2 In frameworks where discretization is part of the regularization of quantum constraint operators, this will 
typically lead to anomalies in the resulting quantum algebra. Hence, one can refer to this phenomenon as classical 
anomalies induced by discretization or discretization anomalies. 
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dynamics defined by the action. This means in particular to allow for a discrete time evolution. 
(A continuous time evolution might be recovered as a symmetry, namely the translation of 
vertices in time direction.) For a discretization based on triangulations such a dynamics has 
been proposed in [H[T0] based on ideas from the consistent discretization program |14j and on the 
so-called tent moves [18], which implement discrete time evolution for triangulated manifolds. 
As a first step we extend these ideas to obtain a canonical description of Regge dynamics. 
However, we cannot expect to find exact constraints in the full theory, as the symmetries are 
broken for curved solutions. Nevertheless, the existence of symmetries for flat solutions should 
have certain repercussions. 

In order to explore these issues, we can consider an expansion of the action on a flat back- 
ground. In this way the calculations become analytically tractable - the full equations of motions 
are so far only solvable by numerical methods. 

The linearized theory, i.e. the theory defined by the expansion of the action to quadratic 
order, has exact symmetries resulting from the null modes of the Hessian on a flat background. 
Hence, the canonical framework for the linearized theory should have exact constraints. We 
will derive these constraints explicitly for an arbitrary 3d triangulation (embeddable into 4d 
flat space), representing the initial data hypersurface. One can show that these constraints are 
Abelian, therefore anomaly free and preserved under time evolution as defined by the tent moves. 
The constraints involve the background geometry in quite a complicated way. Presumably these 
would be hard to rederive by discretizing the continuum constraints directly. Nevertheless, these 
constraints are Abelian, which again follows from deriving them directly from the action but is 
not immediately straightforward to see on the canonical level. Moreover, deriving the canonical 
framework directly from the discretized action gives the momenta an immediate geometrical 
meaning in terms of the discrete geometry. 

For the higher order dynamics one would expect that the symmetries of the action are broken, 
as this is the case for the full dynamics. Indeed, the equations of motion expanded to the lowest 
non-linear order result in consistency conditions on the background gauge, which in the case of 
Regge calculus is associated to the positions of the inner vertices in the flat background solution. 
In other words, a consistent expansion of the solutions (analytically in the expansion parameter) 
is only possible for specific choices of these background gauge parameters. 

The consistency conditions on the background gauge can be rewritten as the condition that 
the quadratic order of Hamilton's principal function (i.e. the quadratic action evaluated on the 
solutions of the linearized theory) has a vanishing derivative with respect to the background 
gauge parameters. In this sense the discretization is fixed such that the dynamics depends 
minimally on the details of the discretization, in this case the choice of background gauge 
parameters. 

In the canonical framework it will turn out that the quadratic order of the constraints 
depends on the background gauge parameter. These could be interpreted as background lapse 
and shift and hence we encounter (lapse and shift dependent) pseudo constraints rather than 
exact constraints. The requirement that the constraints should be preserved under the discrete 
time evolution leads to the same condition on the background gauge parameter as in the covariant 
formalism. 

The structure of the paper is as follows: In section [2] we will give a short introduction to 
Regge calculus and its equations of motion. Subsequently, in section [3l we will discuss the tent 
moves and the associated canonical formalism. In section [5] we explain the origin of the gauge 
symmetries for flat solutions and the relation to the Bianchi identities. The symmetries of the 
flat solution imply that the Hessian of the action evaluated on these solutions will have null 
vectors, which we will examine in section These considerations are essential in order to derive 
the constraints for the linearized theory in section [6J In section [7] and section [8] we will detail the 
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constraints for four- and higher valent vertices. We show in section that the constraints are 
Abelian and, after performing a split into linearized observables and gauge variables in section 
\10\ we consider the dynamics of the linearized observables - that is gravitons - as generated by 
the tent moves in section [TT1 In particular, we will show that the constraints are consistent, that 
is automatically preserved under time evolution. Finally, we discuss the higher order dynamics 
in section [12j We will close with a summary and outlook. 

2 Summary of Regge calculus 

Regge calculus [1J is usually considered on a fixed triangulation T of a space-time manifold. We 
will consider 4d triangulations which are built from 4-simplices. A 4-simplex has five tetrahedra 
r, ten triangles A, ten edges e and five vertices v as subsimplices. Two 4-simplices are glued 
together by identifying a tetrahedron from each 4-simplex with each other. Thus, a tetrahedron 
(in the bulk of the triangulation) is always shared by two 4-simplices, whereas, for instance, a 
triangle can be shared by any number (higher than two) of 4-simplices. 

The variables appearing in the Regge action, which defines the equations of motion for Regge 
calculus, are usually given by the edge lengths {l e } e& j-, for other choices see [191 [20] . These 
variables completely specify the (piecewise linear) geometry of the triangulation. In particular, 
from the edge lengths one can compute any 4d dihedral angle 6^, which give the (inner) angle 
in the 4-simplex a between the two tetrahedra sharing the triangle A. These dihedral angles in 
turn determine the curvature of the triangulation: 

Consider a triangle A in a 4d triangulation. This triangle is shared by several 4-simplices. A 
(Levi-Civita) parallel transport of a vector from one 4-simplex to the next around the triangle 
results in a rotation of this vector by the so-called deficit angle ea = 2-7T — X^ooA @A i n the 
plane perpendicular to this triangle. The deficit angle measures the curvature concentrated at 
the triangle. 

Accordingly, the (Euclidean) Regge action, as a discretization of the Einstein-Hilbert actiorH 
Seh = ~\ J y/gRd 4 x, is given by 



where A a denotes the area of a triangle A. If there is a non- vanishing boundary the boundary 
term 



has to be added to the action in order to make the boundary value problem (with prescribed 
edge lengths on the boundary) well defined. Here ipA = kn — X^coA ^a ^ s the extrinsic curvature 
angle. The value k is determined by how many pieces are glued together at the triangle A in 
question. Usually only two pieces are added, in which case k = 1 to ensure that the actions 
for the two pieces add up correctly to the action for the glued triangulation. Sometimes more 
pieces are glued together, in this case we will use k = for the additional pieces. 

To obtain the equations of motion one has to vary the Regge action (|2.ip with respect to 
the length variables. Here the Schlafli identity is instrumental. The Schlafii identity 




(2.1) 



Acbulk 




(2.2) 



Acbdry 




(2.3) 



Acit 



relates the variatons 59^ of the dihedral angles in a 4-simple3^|. 



3 We work in units with c = 8nG Newton = 1. 

4 This identity can be generalized to any n-dimensional simplex. 
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The Schlafli identity can be understood to be analogous to the result that for the variation 
of the Einstein-Hilbert action — J ^fgg ah R^d^x the term with the variation of the Ricci tensor 
leads to a total divergence. Indeed, also for the Regge action the variations of the deficit angles 
lead to a contribution only from the boundary, which are annihilated by the variation of the 
boundary term (12.2p . The resulting equation of motion obtained by varying the length of the 
edge e in the bulk is 

• (2.4, 

ADe 

A special kind of solutions are flat triangulations, for which all deficit angles vanish, eA = 0. 
Flat solutions, however, can only appear for specific choices of the boundary lengths (if the 
triangulation of the boundary is sufficiently complicatecH) . 

Note that if the boundary lengths allow for a flat solution and the triangulation contains 
vertices in the bulk, this solution is not unique. That is, other flat solutions can be produced 
by translating the inner vertices in the embedding 4d flat space and changing the lengths of 
the inner edges adjacent to these vertices accordingly. On the other hand, one expects that 
for boundary data inducing curvature, the solutions are unique. A family of examples was 
considered numerically in [TO] and uniqueness of the solutions was found. 



3 Canonical formalism 

As mentioned, the edge lengths are not uniquely determined for boundary conditions such 
that the Regge equations admit flat solutions. There is rather a gauge freedom, which can be 
understood from the choice of the exact position of the inner vertices in the flat triangulation. 
Such a gauge freedom is usually accompanied by constraints in the canonical formulation of the 
theory. Here, however, we encounter the situation that only a particular set of solutions, namely 
the flat ones, exhibits an exact form of gauge freedom. We will show in section [6] that this leads 
to constraints in the canonical formulation of linearized theories around flat solutions. Later on, 
we will discuss the repercussions for the full dynamics of the canonical theory. 

To define the canonical formulation along the lines of [6j [TO] we will employ the so-called 
tent moves, introduced in [18j . These admit the advantage to allow for a local evolution of the 
hypersurface on which the canonical data are defined, without changing the connectivity of its 
triangulation. Therefore the number of edges and hence variables remains constant under the 
discrete time evolution. Note that although we use the tent moves to derive the (linearized) 
constraints it will turn out that these are independent from this construction. The constraints 
are conditions on the canonical data, so that these can be consistently evolved. Alternatively, 
(by evolving backwards) one can see the constraints as describing canonical data that can arise 
by an evolution leading to the hypersurface in question. 

To define a tent move consider a 3d triangulation S n , which can be thought of as a trian- 
gulated Cauchy hypersurface with time label n. We will assume that this Cauchy hypersurface 
is a (piece of a) boundary of a 4d triangulation, whose inner edge lengths satisfy the Regge 
equations. Pick a vertex v n in the Cauchy surface and define a new vertex u n +i, which will lie 
in the evolved Cauchy hypersurface S n+ i. Connect both vertices with an edge, which will be 
called 'tent pole'. Denote all other vertices in S n which v n is connected to by 1, . . . , N . Connect 
also v n+ \ to the 1, . . . , N by edges. Furthermore, we will have a tetrahedron T(v n+ iijk) (with 
vertices v n+ i, i, j, k) in the evolved hypersurface for every tetrahedron r(v n ijk) in E n . 

5 There are special (simple) types of boundary triangulations, for instance the boundary of a 4-simplex, for 
which flat solutions are generically possible, i.e. for generic choices of the boundary lengths. The reason is that 
these 3d triangulations can always be embedded into 4d flat space. 
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Hence, the triangulations of the two Cauchy surfaces are the same. The analogous evolution of 
a tent move in 3d is depicted in figures [1] and 



V n +1 




Figure 1: The tent move in 3d applied to a vertex v n in a 2d Cauchy hypersurface. 

The evolution can be thought of as gluing a certain piece of 4-dimensional triangulation onto 
the hypersurface. This 4-dimensional triangulation 7~ n consists of 4-simplices a(v n v n +iijk) for 
every tetrahedron T(v n ijk) in S n . Note that th.6 tent move only involves stctr( , u ri .), the (3d) st&r 
of v n in Yl n . The star of a vertex is the union of all simplices having v as a subsimplex. 

Through this gluing-on we have obtained an additional (N + 1) inner edges, namely the 
edges e = e(v n i),i = 1, . . . , N and the tent pole t = e{v n v n+ \). We will rewrite the equations of 
motion for these edges into a canonical form. To this end denote by 

S n = - Y, [27T- Y, 9 a] 



AcT„ 



- E a a u-y°a - Y a a u-y^ 

ACstar(i! n ) ACstar(u n +i) 

Y ^a(-E^) (3-1) 

the Regge action for the added piece of 4d triangulation T n (with boundary terms). With 

o o 

A G star(u n ) (or A G star(u n+ i)) we mean triangles that are in star( , y n ) but are not part of 
star(v n+ i) (or vice versa). There are also triangles which are part of both S n and S n +l- If one 
performs several consecutive tent moves at the vertices v n ,v n +i, . . . , then these triangles are 
part of each of the triangulations 7~ n ,7~ n +i, ■ ■ ■ ■ Hence, we choose the associated boundary term 
without any factor of tt, as we cannot say how many pieces T are added. (Also if tent moves 
at neighboring vertices are performed then the action associated to these moves provides the 
necessary factors of tt for these triangles.) 

With S n ~x we will denote the action (again with boundary terms) of the original 4d triangu- 
lation without the piece T n . (Alternatively, one can assume that a tent move at v n -\ has already 
been performed. Then S n -\ is the action associated to T n -\. Again, this does not matter for 
the equations of motion.) The equations of motion can be written as 

dS n 

9t n 

dS n —\ dS n , . 
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where by t n we denote the length of the tent pole t = e(v n v n+ i) and l n is the length of the edge 
e = e(v n i),i = 1, . . . , N. Using S n as a generating function, we define the momenta canonically 
conjugate to l e n , t n , t n+1 by 



Pt 



Pt 



n+l ._ 



_dSn 

dt n 
dt n +i 



Pe 



Pe 



dSn 
~ 81* 

dl 



(3.3) 



n+l 



Note that the momentum p™ +1 identically vanishes as S n does not depend on t n+ i. The equations 
of motion (|3.2[) are now simply given by 



Pi 
Pi 



dt n 
dS n -i 



dS.. 



dt 

BS n 



dl n dl* 

and thus reproduce the Regge equations of motion (13. 2p . 



-=p? = o 

Pe , 




(3.4) 



Figure 2: The local tent move evolution of a vertex v n in 3d. 

Apart from the edges e adjacent to v n there are more edges b in (the boundary of) the 3d 
star of v n . The lengths of these edges do not change under a tent move at v n , however, if one 
performs tent moves at neighboring vertices one has to transform the momenta associated to 
these edges. This transformation is dictated by the requirement to reproduce the Regge equation 
of motion for all tent moves. This is achieved by defining a generating function (of mixed type) 

F(Jni ^n+li^ni ^n>Pb + ) = ~ ^~]^n Pb~*~ + Sn(J>ni ^n+li *m O • (3-5) 



The transformations for the variables associated to t,e do not change from f|3. 3|) . For the edges 
b in the boundary of star(u n ) we obtain 

l b - dF -l» 

'n+l - g p n + l ~ 

Pb - ~W n ~ Pb 8% • (3 - 6) 

As we have l n = l n+ i we will often use just l b for these variables. 

The canonical transformation (|3.3l I3.6|) (it is canonical because we defined it via a generating 
function) defines the evolution of the canonical data from time step n to time step (n + l). Here 
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the equation of motion pf = is analogous to the primary constraints appearing in continuum 
(canonical) general relativity pj^ = 0, which imply that the momentum conjugate to the lapse 
function vanishes. The momenta at time step n as defined by (|3.3|) are explicitly 

N aA 



it Ql fc A(i) n v n+ ii) 



i=l 



Pe = 2^ ~di^ lJjA + — d~F n eA (^+^) (3 - 7) 

Agstar(ii n ) 

with e = e(v n i) - the edge connecting v n and the vertex i. For the new momenta we obtain 

n + 1 ^ 8A A dA^ v n v n + l^ 

Pe ~ ~ ~Q^e ^ A Qje e A(v n v n+1 i) 

o n+1 
A£star(ti n +i) 

n+1 n ST^ 9A A ® Aa I ^ Aa I (n o\ 

p b = p b - 2^ ~dF^ A ~ 1^ ~~dF^ A ~ ~dF^ A - (3 - 8) 

Aestark) Aestar(; n+1 ) Ae S tar( Vn )nstar(„ n+1 ) 



4 Gauge symmetry and contracted Bianchi identities 

In continuum general relativity the contracted Bianchi identities, 

V a G ab = , (4.1) 

for the Einstein tensor G a b = R a b — ^RSab can be derived from the invariance of the Einstein- 
Hilbert action under diffeomorphisms, see for instance |2H I22j . On the other hand, the contracted 
Bianchi identities are geometrical identites, which follow from the properties of the curvature 
tensor. 

From the contracted Bianchi identities it follows that not all ten of Einstein's field equations 
are independent. The identities V a G a b = provide four differential relations between the field 
equations; in other words, the field equations are not fully independent from each other. Given 
sufficient initial data, the evolution of the ten metric components is therefore not completely 
determined. This is related to the freedom of choosing coordinates and therefore ultimately to 
general covariance. 

Similar arguments can be made for Regge calculus [U El El [23] , in which the Bianchi identities 
hold as geometrical identities, giving relations between finite rotation matrices |231I24|. There is, 
however, a difference to the continuum case, namely that the relation to the equations of motion 
can only be made approximately, only valid for small deficit angles (and further assumptions on 
the 'fatness' of the simplices [5]). This means that the equations of motion will be dependent 
only in this approximation. However, this approximation turns into an exact identity for the 
linearized theory on a flat background. 

In the following discussion we will adapt some arguments from [5] to the linearized theory 
and clarify the origin of the degeneracy of the Hessian of the Regge action which we will discuss 
in section [5j In a triangulation of flat space one can displace the vertices in the embedding 
flat space without changing the deficit angles (which are vanishing). The induced infinitesimal 
change of length variables 5l e is described by the vector fields 5l e =Yf (J = 1, . . . , 4) which can 
easily be computed to 

Yf = (4-2) 
V E e ■ E e 
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where Bj is a basis in the embedding flat 4d space and E e are the 4d vectors for the edges e 
adjacent to a given vertex v (all edges are either pointing towards the vertex or pointing away). 
For edges e not adjacent to v the components of Yf are zero. The vector fields Yf act trivially 
on the deficit angles (evaluated on flat backgrounds), that is 

Taking the sum over triangles and multiplying with a factor dA A /dl e ' we obtain 

A 

Note that the Hessian of the Regge action (with respect to the length variables of edges in 
the bulk) evaluated on a flat background is exactly 

d 2 S _ ^ de A 

WdF W di e lflat [ ' 

A 

so that equation (|4.4|) shows that this Hessian is degenerate and that the derivatives appearing 
in (|4.4p commute. Expanding the lengths as l e = (°H e + ey e + 0{e 2 ) we therefore find that the 
equations of motion for the linearized theory are linearly dependent, 

Y il^-Qir-Qi^\fiaty e =° • (4-6) 

A 

This equation is a first order expansion of the 'approximate Bianchi identity' [H [5] 

a ^«0, (4.7) 

expressing that the equations of motion for the full theory are approximately dependent on each 
other. Conversely, equation (|4.7p can be used to show that the Hessian of the Regge action is 
degenerate - equation (|4.4p - when the linearized Bianchi identity holds. 

The fact that the equations of motions are not independent from each other leaves some of 
the variables - four per vertex - undetermined, which explains the gauge freedom appearing in 
the (linearized) theory. 



5 Degeneracies of the Hessian of the action 

The linearized theory is defined by the quadratic expansion of the action around a background 
solution. As the linear terms vanish due to the background satisfying the equations of motions, 
the linearized theory is basically determined by the Hessian of the action, i.e. the matrix of 
second derivatives. Consequently, we need to understand the properties of this Hessian (here 
evaluated on a flat solution) in order to understand the properties of the linearized theory. 

Assume that we have performed two consecutive tent moves so that the triangulation that 
was added in the process is T n -i U T n - This piece of triangulation has one inner vertex v n . 
Consider this piece of triangulation as a boundary value problem, i.e. fix all the edge lengths in 
the boundary S n _i U S n+ i and solve the Regge equations for the {N + 2) inner edge lengths 
{Z^,e = e(v n i),i = 1, ...,N} and t n -i,t n . Furthermore, assume that the boundary lengths 
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admit a flat solution, i.e. a solution for which all the deficit angles at the inner triangles vanish^ 
Now the triangulation T n -i^T n represents a piece of flat space and by moving the vertex around 
inside this piece of flat 4d space and changing the length of the adjacent edges appropriately, 
i.e. according to their embedding into flat space, we will obtain a 4-parameter set of further 
solutions. 

Hence, the extremum of the action corresponding to the original flat solution is not an 
isolated one, rather there are four constant directions at this extremum. Accordingly, and as 
discussed in the previous section, the Hessian of this action with respect to the inner edge 
lengths {1^, e = e(v n i),i = 1, . . . , N} and t n -i,t n has four null vectors Yi, I = 1, . . . , 4, whose 
components we will denote by (Yj"' 1 ,Yj n , {Yf}). 

In the following we will show that from these null vectors we can define null vectors for 
other Hessians that will appear in the dynamics of the linearized theory. To begin with we will 
eliminate the lengths of the tent poles as these function are auxiliary variables. To this end we 
define the 'effective action' 

^(^n-l^n^n+l^ ) := <Sn— 1 (T n -\ (l n _i , l n , I ), l n _i,l n ,l ) + S n (T n (l n ,l n+l ,l ), l, n ,l n+l ,l ) 
= <5n-l(^n-l^n; O + £>n Q n > > O (5-1) 

where we have solved the associated equations of motion for the length of the tent poles, so that 
dS i dS 

-q^ -(T n -l{l n -l, ^n)> ^n-l>O = 0> (^n(^>^n+l)> ^n^n+l) = ^ ■ (5-2) 

The following arguments will show that the {Yf} define the null vectors of the Hessian of 
the 'effective action'. Taking the derivative of these equations (|5.2p with respect to l n we obtain 
the identities 

d 2 S n -\ d 2 S n -i (9T ra _i g 



dt n -idl% dtn-xdtn-i dl% 



d 2 S n d 2 S n dT n 

+ 51 51 o,p = . ( 5 - 3 ) 



dt n dl e n dt n dt n dl e n 



so that 





( d 2 S n ^ 




\dt n _idt n - 


dT n 


( d 2 s n y 




\dt n dt n J 



d 2 S; 



n-1 



dt n - l dl™ 
1 d 2 S n 



dt n dl 



(5.4) 



e 



On the other hand, following from the fact that the Yj are null vectors of the Hessian of 
S = S'n-i + S n (we will apply the Einstein summation convention for the index e) we have 

ye d 2 S tn _, d 2 S . d 2 S 

1 dl*dl* 1 dtn-idl* 1 dtndl* 

ye 9 2 S ^ d 2 S = 

1 dl™dt n ~ l 1 dtn-idtn-i 

dl n dt n dt n dt n 



6 Such boundary conditions can be constructed as follows: Start with the piece of 3d triangulation E n _i and 
embed this into flat 4d space. Now perform the tent move but within the embedding of the flat space. That is, 
after choosing the position of the tip of the tent pole v n +i (and the lengths in E n _i) the lengths of the edges 
from the tip of the tent pole to the vertices i = 1, . . . , N are determined. This gives a boundary that - as it is 
embedded in flat space - admits a flat solution. The boundary problem is solved by just choosing the position of 
the inner vertex v n and connecting this vertex with the vertices in the boundary appropriately. 
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From the last two equations in (|5.5|) we obtain the components Yj" 1 and Yj n as functions of 



Yf 



yt„-i ye 9_Sn-l ( (PSn-1 

1 1 dl^dt^ V^n-l3t n _l 

Yt{X = " y/ ^^(^jt) " (5 ' 6) 

(Here we assume that the second partial derivatives of the action with respect to t n -\ and t n 
do not vanish. This is generically the case for the Regge action.) 

The first equation in (|5.5p together with (|5.4p and (|5.6p can be used to show that the 
components {Yf} constitute null vectors for the Hessian of the effective action S, that is 

ye d*S 

1 dl*dl£ 

d 2 S n -\ ( d 2 S n ~\ \ d 2 S n -\ 8 2 S n ( d 2 S n \ d 2 S n \ 




61*81% dl e n dt n -x \dt n -xdt n -\) dt n ^dl% dl e n dt n \dt n 8t n J 6t n 8l% J 

= / (5.7) 

Later on, in section [TT1 we will need this relation in the form 

d 2 S n -i d 2 S n 

Yl di^M + Yl di^M-° • (5 - 8) 



Similarly, we have 



d 2 S n -i d 2 S n _ 



where l b is the length of an edge which is contained in both three-dimensional star(f n _i) and 
star(v n ). This equation follows if we consider the Hessian associated to a larger boundary 
problem, also including l b as free variable. Since a translation of a vertex v n only affects the 
lengths of the edges adjacent to v n this Hessian still has the null vectors Yj (with the components 
i / e,t n -i,t n vanishing), in particular, 

Y! 4Sr b+Ytrl di^ +Ytin ^3r b = • (5 - 10) 

Together with the equations (j5.6j) for the components Yj™' 1 and Yj n , we obtain the result 
asserted in equation (I5.9p . 

Next, we will show that the null vectors Yf are also left or right null vectors of 

d 2 S n d 2 S n -\ , . 

or 7^ — ^ , (5-11) 



respectively. Assume that we extremize the action S with respect to the variables Z^, fixing 
the variables Z^_ x and Z^ +1 . Calling the corresponding solutions L^(Z^_ 1 , Z^ +1 ), we obtain the 
identities 

6S 
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Differentiating these equations with respect to or l^-i results in 



a 2 s dLi' j d 2 s n 



respectively. This proves the claim. 

Finally, let us note that the components Yf can be computed from the lengths of the edges 
in the 3d star(u n ) only. As was noted in section we have Yf = Bj ■ E e /\E e j\ where Bj 
is a basis in 4d flat space and E e are the 4d vectors for the edges e, which one can obtain 
by embedding star(u n ) into 4d flat space. From a counting of variables argument, which can 
be found in section [8l one can deduce that the lengths in star(t>„) determine this embedding 
uniquely (modulo translations and rotations). Hence, the components Yf are determined by 
and l b . As we will see in the next sections, the same holds for the constraints of the linearized 
theory, i.e. these depend only on the (fluctuation and background) variables associated to a 
particular (discrete) time. 



6 The linearized theory 

Here we have to consider the tent move equations ()3. 3f) to linear order. That is, we expand all 
lengths l e of the triangulation around a flat solution in a small parameter e and keep only 
terms linear in e in the equations. Later on, we will also discuss an expansion to higher order 
in s. 

Among the possible variations of the lengths are also the ones in flat directions. These 
variations are solutions of the linearized equations of motion. As solutions to linear equations 
of motions are additive these variations can also be added to solutions representing linearized 
curvature excitations without changing the boundary data. All the solutions of the linearized 
theory will therefore exhibit gauge symmetries. In general, the linearized theory will inherit the 
gauge freedom of the background solution. 

In a tent move we encounter three different kind of edges. Firstly, there are edges b in 
the intersection of the two Cauchy surfaces £„ n defined by the tent move. These edges 

are not dynamical (for this specific tent move) and we will call the associated length variables 
l b = (o)/ 6 _|_ £ y b _|_ 0(e 2 ). Secondly, there are the 'dynamical edges' in the Cauchy surfaces S n 
and S n +i adjacent to the evolving vertices v n ,v n+ i. The associated length variables will be 
l e = (°)/ e + ey e + 0{e 2 ). Finally, there is the tent pole - the only bulk edge - with length 
t n = {0) t n + ex n + 0{e 2 ). 

Furthermore, we expand the momenta conjugate to 1^ and t n , that is = +£71"™ + 0(e 2 ) 
and p n = (°'Pt + £7r™ + 0(e 2 ). To linear order in e we obtain from the first two equations in 
(|3.3p . which define the momenta and p™, 

n _ d 2 S n e i d 2 S n e / d 2 S n d 2 S n b . . 

^ " ~mM Vn ~ di- n dii +l Vn+l ~ m^dT n Xn ~ dW bV ' 

n _ d 2 S n e d 2 S n e d 2 S n d 2 S n ^ 

nt ~ ~dt^dF n Vn ~ dt n di e n+1 Vn+1 ~ £ dt^df n Xn ~ dt^ y 

= , (6.2) 

where it is understood that the derivatives of the action S n are evaluated on the flat background 
solution. 
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Solving equation (|6.2p for the tent pole variable x n and using this in (|6.1|) . we obtain 

< = -M: e ,yi +l -N: e ,yi-N: b y b , (6.3) 



where 



N n , 

x, ee' 



a 2 s n d 2 s n ( d 2 s n Y 1 d 2 s n d 2 s n 



dl'dl^ dl n dt n \dt n dtj 8t n dl* +1 8l*8l* +l 



d 2 S n 


d 2 S n 


/ d 2 S n ^ 


f 1 d 2 S n 


d 2 S n 


81*81* 


8l n dt n 


\8t n 8t n/ 


I 8t n 8l* 


81*81* 


d 2 S n 


d 2 S n , 


' d 2 S n \ 


1 d 2 S n 


d 2 S n 


dl e n dl b 


8l*8t n ' 


ydtndtnj 


dt n dl b 


dl e n dl b 



eb di n di b di n dt n \dt n dt n J dt n di b di*di b ' y ' J 

From the discussion in section [5] we know that the first matrix M"/ has four (left) null 
eigenvectors Yf,I = 1,...,4. Contracting equation (|6.3|) with the left null eigenvectors we 
find four relations involving only momenta and configuration variables at time step n, that is 
constraints 

Cj = Yf^ + YfN2 e ,y* + YfN2 b y b = . (6.5) 

These constraints have to be satisfied by the canonical data of the linearized theory in order to 
give rise to a solutions of the equations of motion. Although the constraints involve only the 
momenta and configuration variables of linear order, variables describing the flat background 
solution might appear, which also involve data at the time step (n + 1). However, these could 
be replaced-using the equation of motion for the background-by the configuration variables and 
momenta at time step n as well as background lapse and shift (describing the gauge freedom of 
the background solution). But, as we will see, this issue will not arise. That is, the constraints 
will only involve the background variables at time step n and not depend on the position of the 
vertex v n +\ of the flat background solution. 



7 The constraints at a four— valent vertex 

Here we will discuss the constraints ()6.5p in more detail and derive a more explicit expression. 
We will start with constraints derived from the tent move at a four-valent vertex. A four-valent 
vertex in the 3d boundary of a 4d triangulation can be identified as a vertex of a 4-simplex of 
this 4d triangulation. For a four-valent tent move one carjj| construct a flat solution]! Indeed, 
having given the four edge lengths £™ and e = e(vl), . . . , e(v4) in addition to the six edge 
lengths If,, b = e(12), . . . , e(34) we can construct a solution by taking two 4-simplices <j(v n 1234) 
and cr(v n+ il234) with the appropriate edge lengths and gluing these together along r(1234). 
Connecting v n with v n +i, we will obtain the length of the tent pole t n . All the deficit angles at 
the triangles hinging at the tent pole vanish, hence the Regge equation associated to the tent 
pole is satisfied. The analogous construction in 3d is depicted in figure [3) 
With the solution T n (l n ,l n+1 ,l b ) of the tent pole equation, 

t 8S n s-^ dA A 

ADt n 



7 Under the condition that the triangle inequalities are satisfied in the construction below . 

8 There might also exist exceptional cases where solutions with curvature are possible, see the discussion in 
[10] . These seem to be discretization artifacts, however, and we will ignore these kind of solutions. Moreover, also 
the flat solutions might be ambiguous, in particular there will be future directed and past directed solutions. We 
will always choose the future directed solution. 
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Figure 3: Construction of the 3d three-valent analogue of the 4d four-valent tent move where both 
tetrahedra are oriented in the same (future) direction. 

we can define the following functions of in^n+i an d l h 

n n — <n n , — 

Pb '■= Pb\t=T n (lZ,,ie,,,l<>) := oT t • ( 7 - 2 ) 

Taking the derivative of equation (|7.ip with respect to , we obtain 

dT n= f d^n y'd^n 

dll \dt n dt n ) dtj* [ ' ' 

and can conclude that the matrices N™, e , N™, b defined in ()6.4p can be expressed as 

dit ee ' dii eb ■ { 1 

As discussed above, for the four-valent tent move the deficit angles at the triangles hinging 
at the tent pole vanish, hence the expressions for pf l ,p b simplify to 



Ac£„ 

di h 



"" E -^M^CiA ■ (7-5) 



AcS n US n+ i 

Here ^A^™'' ' ^ are ^ ne extrinsic curvature angles for the boundary triangles of a piece 
of flat triangulation. For the situation depicted in figure these angles are given by 

V>A = -7r + A (C'/) forAci 
ipA = n-OA(l e n+i,l b ) forAcE n+1 

^A = -O A (ln +1 ,l b )+0A(lnJ b ) for A C £„ D £ n+ i , (7.6) 
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where #a(0 is the dihedral angle on a 4-simplex a with edge lengths I. 

In order to compute the matrices iVp e and N™, b , note that taking the derivative with respect 
to l!^ will annihilate the dihedral angles #a(^+d l b ) m (|7.7p . so that everything can be expressed 
with dihedral angles of the simplex a with edge lengths (7^, l b ): 

^ = ^E^(--^'/)) 

n Ac<7 n 

= Jr£^-M^ 6 )), (7.7) 

where we used the fact that the right-hand sides of both sides is a double derivative of S a = 
Y2a ^-a(tt — #a)- Formulae for the derivatives of dihedral angles can be found in [9]. 

The same applies to the computation of iVp b , moreover, note that d/dl^A^j^ = as 
^A(ijfc) om y depends on the lengths l b . Thus, we may write 

n \ACa,A=A(v ij) ACa,A=A(ijk) 
ACcr ™ 

Note that the dependence of the matrices N^, e ,N^ b on the lengths drops out completely. 
That is, we can compute the constraints from the (background) geometrical data on S n only. 
To finally compute the constraints 

d = Yf{i% + iV™,/ + N eh y b ) (7.9) 

we need to determine the four vector fields Yf . As we have four vector fields and four edges 
e = e{vQi) we can identify the indices I and e and define Y5 to be the vector field that translates 
the vertex vo orthogonal to the other three edges e ^ e'. It is easy to see that (assuming 
normalization) we have Y5 = <5f/. This determines the first order constraints for a four-valent 
vertex to 

c < = < + w £ 9 ~W (7r - 9 ^ lb ))y e n + w E ^ (* - Mtf . W • (7-io) 

n ACcr n ACcr n 

Note that these constraints coincide with the first order of the full constraints derived in 

ci un =p n e + Y, d 4^^- e ^ 1 "^ ■ (7 - n) 

ACff n 

These constraints for the four valent vertex can be generalized to the 'flat sector' of 4d Regge 
calculus, which are a special class of triangulated hypersurfaces that only allow for evolution 
leading to flat 4d space, see [T7] . 



8 The constraints at higher valent vertices 

Next, we will derive the constraints at higher valent vertices. The discussion will be in many 
aspects parallel to the one in the last section [71 There is one important difference however, 
which is that solving the tent pole equation for higher valent vertices can also lead to solutions 
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with non-vanishing deficit angles. Nevertheless, we will see that again for the computation of 
the constraints we will only need the geometrical data of the hypersurface 

To this end, note that a tent move at a vertex v only involves the 3d star of this vertex, 
i.e. all simplices (with their subsimplices) that share the vertex v. As we perform an expansion 
around flat space the configuration we are considering must be embeddable into flat 4d space. 
Indeed, there are as many edge lengths in the 3d star of a vertex as one needs to determine 
an embedding into flat 4d space modulo translations and rotations [8]. Firstly, we will count 
the number of edges in the star of an iV-valent vertex. In addition to the N edges adjacent to 
v we have E edges in the boundary of the star. The piecewise linear manifold condition [25] 
ensures that this boundary is topologically a 2-sphere. For the number of edges E, the number 
of triangles T and the number of vertices V in a triangulated 2-sphere there are two relations, 
the Euler theorem T — E + V = 2 and the relation 3T = IE. Hence, the number of edges in the 
2-sphere is E = 3^ — 6. The number of vertices is V = N so that the overall number of edges in 
the 3d star is 4iV — 6. On the other hand, if we embed the 3d star of the iV-valent vertex v into 
4d flat space we have to choose AN coordinates for the N vertices. Modulo the six 4d rotations 
this will also give AN — 6 parameters. Hence, one can expect that the lengths of the edges in 
the 3d star uniquely determine an embedding into flat 4d space. (There might occur discrete 
ambiguities, however, these are fixed by the flat background solution under consideration.) 

To derive the constraints we will again use the momenta introduced in (|7.2p in order 

to compute the matrices and N™, b . These momenta involve the solution of the tent pole 
equation (|7.ip which, however, can now lead to solutions with non-vanishing deficit angles. That 
is, we have 

Pn\flat 
V \flat 

only on data l^ +1 ,l b ) which lead to a solution of the tent pole equation with vanishing deficit 
angles. (Fixing l b ), that is the geometry of the 3d star S n , there is generically a 4-parameter 
set of lengths l^+i sucn that one can obtain a flat solution.) 

This is, however, sufficient to compute the contraction of the matrices A^p , iVI). with the 
vectors Yf, which is all we need to determine the constraints. According to (|7.4|) . this contraction 
corresponds to the derivatives of the momenta p in the direction of Yf: 

Y?^jjr = -YfNX , Yf 'M = - Y f ' N% . (8.2) 

As explained in section HJ these vectors correspond to translations of the vertex v n and the 
induced change of lengths such that the triangulation remains flat. Hence, we can still use 
the expression (]8.ip to determine the derivatives in 'flat directions'. As for the 4-valent vertex 
we have again 

V>a = -7T + e A (l e n ,l b ) forAci 
^A=ir-0 A (l e n+1 ,l b ) forAcIWi 

^ A = -e A (i e n+1 ,i b ) + e A (i e n ,i b ) forAcs n ns n+1 , (8.3) 

where we now generalized the angles 6 A to the dihedral angle between the two tetrahedra 
sharing the triangle A of the 3d star of v n or i> n +i, respectively, embedded into 4d flat space. 



AcE„ 

E ^Ml e n,l G n + l,l b ) , (8-1) 
AcS„US„ + i 
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As discussed above, these embeddings, and hence the dihedral angles, are determined by the 
edge lengths or respectively. Again, the dihedral angles 9/^(l^ +1 ,l b ) drop out 

after taking the derivatives with respect to . As explained in section [5j the vectors Yf can be 
determined as functions of 1^, l b ■ The linearized constraints take the form 

0A(l e n ,l b ))y e n 

A (l*,l b ))y b . (8.4) 

Also here the linearized constraints do not depend on the lengths in the background solution 
at the next time step, in particular they do not depend on background lapse and shift. In this 
sense the linearized constraints are independent from the tent move construction: they refer 
only to the data on the Cauchy surface in question, hence we do not need to specify any tent 
moves in order to state the constraints. 

The constraints (|8.4p also generate the expected gauge transformations. The constraint Cj 
is expected to generate the change of coordinates induced by translating the vertex v n in the 
direction determined by Yf (whose components describe the induced change of length variables). 
Indeed, 

{y e n ,d} = Yf . (8.5) 

To determine the infinitessimal change of the momenta, remember that translating a vertex 
in the embedding flat space does not change the flatness of the configurations. On such flat 
configurations the momenta are given by (|8.1j) . (For the boundary edges with index b formula 
(|8.ip gives the part that does depend on the length l e n of the edges adjacent to v n .) Therefore, 

{K^Ci} 
{vr b ,C/} 

reproduces the correct transformation behaviour for the linearized momenta. Hence, requiring 
that the constraints Cj generate the change of variables induced by vertex translations in the 
direction Yf gives an alternative derivation of the formula (|8.4p . 

8.1 Example: five— valent symmetry— reduced vertex 

As an example we consider a tent move at a five-valent vertex. To simplify the situation, we will 
consider a 'symmetry-reduced' set-up, also used in jlOj , so that we only have two dynamical 
length variables, a n and b n , to deal with at each time step. 

The geometry of the 3d star(-u) is illustrated in figure 2] and given as follows: 
As v is five-valent we have five further vertices which we will denote by 1, . . . , 5. We will assume 
that we have six tetrahedra with vertices 

«124, fl34, v234, vl25, «135, v235 . (8.7) 

Accordingly, we will have nine triangles of the form A(vij) with i,j = 1, . . . 5 in this trian- 
gulation, five edges of the form e(vi) and nine edges of the form e(ij) (all possible ordered 
combinations of i,j G {1, ... 5} with the exception 45). 



r ve'n, V e' ° \" 9A A 

AcE„ 

, V ef 9 sr 9Aa U 

n AcS n 



- Y r'A E i^(-^A(C/ fe )) 



AcE„ 



(8.6) 
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Figure 4: Illustration of the symmetry-reduced 3d star(u) of a five-valent vertex v, consisting of the six 
tetrahedrar(ul24), r(ul34), t(v234), t(«125), r(vl35) and r{v235). 

The symmetry assumption includes to set all the lengths of the boundary edges e(ij) to 
1 l!l (and we will neglect variations y b at these boundary edges) and setting equal l e ( m > = a, 
i = 1, 2, 3 as well as = = b. 

The 4-simplices involved in the tent move are then all of the same type a{voViijK) where 
vo, v\ denote the two vertices of the tent pole (at time steps n = 0, n = 1, respectively), i, j take 
values in 1,2,3 and k, in 4, 5. We will denote by 

6>°, A° a the dihedral angle and the area of the triangle A(voij), 

8 b , A® the dihedral angle and the area of the triangle A(voin), 

01, the dihedral angle and the area of the triangle A{voV\i), 

8 b , A\ the dihedral angle and the area of the triangle A(vqV\k\ 

d\, A\ the dihedral angle and the area of the triangle A(v\ij), 

8 b , A\ the dihedral angle and the area of the triangle A(v\in) 

8, A the dihedral angle and the area of the triangle A(zj'k), respectively . 

The Regge action for one time step is given by 

S = -3At(2n-49-)-2A b t (27r-39 b t )-3A° a (7r-29 a )-6A b (7r-28 b ) 

-3A 1 a (n-2e° a )-6Al^-29 b )-6A(-8) . (8.8) 

We define the momenta to be 

o dSo dSo , , 

P « = -^' Pb = ~Wo ■ (8 ' 9) 

On flat configurations we have 

where for the exterior angles we can write 

^° a :=ir- 29 a = -ir + 2d a (a ,b ), ^ = tt - 28° b = -tt + 9 b (a , b Q ) . (8.11) 

Here 8 a (a, b) and 8b(a, b) are the dihedral angles in a simplex a{vijkn) with e(ij) = 1, e(vi) = a 
and e{vn) = b at the triangles A(vij) and A{vin), respectively. The reason for these relations is 



9 This is possible since the vacuum Regge equations are invariant under global rescalings. 
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that the 3d star(u ) under consideration can be constructed by gluing two 4-simplices <t(uo1234) 
and <r(fol235) together along their common tetrahedron t(vq123). The 3d star(t>o) is given 
by all the tetrahedra except for r(1234), r(1235) and r(t>ol23). If we glue the triangulation 
corresponding to the tent move (with six 4-simplices) and these two 4-simplices along the 3d 
star(fo) together we can use flatness of the deficit angles at the triangles A(voij) and A(voin) 
to conclude (|53D|) . 

Hence, the momenta for flat configurations can be expressed as functions of the configuration 
variables ao,feo only. Similarly, the components of the vector Y a ,Y b induced by displacing the 
vertex vq in the embedding flat spacetime can be expressed as functions of the configuration 
variables ao, &o only. Note that due to our symmetry requirements there is just one such vector. 
A displacement of the vertex vq can only change the data associated to the star of vq. For the 
complex of the two glued simplices it should not change the exterior angle V'A(i23) a ^ the triangle 
A(123) between the tetrahedra r(1234) and r(1235). This gives one relation between the two 
components of the vector Y which can be computed to 

y a = 4__I, y b= 4 + bl-l (812) 

a 26 

and satisfies 

^ a ^A(i23) + Y b —^ A(123) = . (8.13) 

The linear constraint is according to (|8.4|) the projection of the expression (18.10P for the flat 
momenta onto the vector Y, that is 

C = Y a ir° n +Y b ir° h + 



f) 8 ( f)A° BA° 



8A° 

6-^(K-e b (a ,b )))y b o ■) cs-ii) 



9 Constraint algebra 

Gauge symmetries lead to constraints which in turn generate the gauge transformations in the 
canonical framework. The fact that the gauge transformations form a group is reflected in the 
first class property of the constraints, i.e. the Poisson bracket of two constraints should give a 
combination of constraints or vanish. 

In this section we will show that the linearized constraints (|8.4p are indeed Abelian and 
therefore first class. (Constraints that are linear and first class have to be Abelian.) As we will 
also consider constraints based at different vertices v, v' we will change the notation slightly, in 
particular we will omit the index for the time step n and introduce an index v for the constraints 
v Cj and the vectors v Yf based at the vertex v. We define v Yf = if e is not adjacent to v. 

We are considering an expansion around flat 4d space, hence we can assume that the trian- 
gulated hypersurface S (with the background edge lengths l e ) is embedded into flat 4d space. 
Thus, there is some suitable flat 4d triangulation T such that £ is (part of) the boundary B{T) 
of 7". Furthermore, we will use the index e for any edge in S, i.e. not only for the edges adjacent 
to v or v' , and therefore not use the index b anymore. By star(f) we denote the 3d star of the 
vertex v in S. Writing all summations over the index e explicitly, the constraints fj8.4f) are now 

^ = 5> /7 r e+ £ E'YfA £ ^(vr-W . (9.1) 

eDv e'Cstar(i)) eDn ACstar(i;) 
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The Poisson bracket between two constraints is given by 



'f'CeCstar(ti) e'Dv ACstar(u) 

- E E S 8 -w { '-^ ■ (9 - 2 » 

uCeCstar(u') e'Dv' ACstar(u') 

We will show that the two terms on the right hand side of (|9.2p cancel each other and hence 
the constraints commute. To this end we will prove that both terms are second derivatives of 

S T := A A (ir-8 A )+ ^ j4a^a (9.3) 

Ac£>(T) AcX(T) 

contracted with '"Yf and v Yj. Here A C B(T),I(T) denote triangles in the boundary and bulk 
of the 4d triangulation T, respectively. Using the Schlafli identity the first derivative evaluated 
on a flat configuration amounts to 

AcB(T) 

As the second derivative is contracted with a vector '"Yf or v 'Yj along which the configuration 
stays flat, we can still use the expression for the first derivative. Finally, note that 

e e' ^ ' 

is only non-vanishing, if the triangle A is both in star(u) and star(w'). Firstly, the derivative of 
the area A A is only non-zero for e C A, for the second derivative we additionally need e' C A. 
Secondly, the derivative of the dihedral angle #a is contracted with a vector that arises from 
displacing the vertex v' in the flat 4d embedding space. Under such a displacement only dihedral 
angles associated to triangles in the star of v' are affected (as only edges adjacent to v' change, 
consequently, only the normals to the tetrahedra in the star of v' , and the normals determine the 
dihedral angles). This shows that the second derivatives of Sj- contracted with ~"Yf and v Yf in 
the two possible ways give, indeed, the two terms in the Poisson bracket (j9.2|) . which therefore 
vanishes. 



10 Observables 

We have seen that in the canonical framework for linearized Regge calculus constraints appear. 
These have an Abelian algebra and, moreover, generate gauge transformations that correspond to 
translating the vertices of the hypersurface in the embedding flat (background) space time. Thus, 
not all variables are physical observables, that is invariant under gauge transformations. Gauge 
invariant variables are those that do not change under vertex translations and are associated to 
(linearized) curvature excitations. 

Considering the dynamical variables for a tent move at some A^-valent vertex v (we again 
omit the time step index n) we define a canonical transformation 

/ = (T- 1 f e y e 

vr r = T r e vr e , (10.1) 
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where we divide the index set T = 1,...,N into two sets / = 1,...,4 and a = 5,...,N, 
corresponding to gauge variables and gauge invariant variables, respectively. For the y a and 7r a 
to commute with the constraints we need that 

Y?(T- 1 )° = 0, YfN ee ,T^ = . (10.2) 

Choosing Tf = Yf ensures the first equation in (|l(J.2p . moreover, it then follows that the 
7Tj = Yfir e coincide with the momentum parts of the constraints. In particular, with the second 
condition in f)10.2|) 

d = it j + YfN ee iYp y 1 ' + YfN^j, y b , (10.3) 

that is the constraints only involve the gauge variables y 1 ,ttj. 

Partial and complete observables [26] are a general tool to compute physical observables and 
can also be applied here. Select four linearly independent edges (with respect to the geometry 
of the background solution) among the N edges adjacent to v. We will split accordingly the 
index set e = 1, . . . , N into E = 1, . . . , 4 and P = 5, . . . N. The Poisson brackets between y E 
and the constraints give then an invertible matrix 

{y E ,C I } = Y I E . (10.4) 

(Note that the vectors Yf could have been chosen in a way such that Yf = Sf , see the expla- 
nation in section O) To any linear phase space function / we can then associate an observable 
Ff defined by 

F f = f-{f,C I }(Y- 1 ) I E y E . (10.5) 

We can similarly proceed with any set of four momenta tt e ', such that the matrix {ir E ,Cj} is 
invertible. A geometrical interpretation of these observables will appear in |27j . 

In the following we will discuss the structure of the phase space for the dynamical variables 
of a tent move at an iV-valent vertex v. As we have four constraints the constraint hyper- 
surface is (2iV — 4)-dimensional. An A^dimensional submanifold of this hypersurface is given 
by configurations (here linearized length and momentum variables) leading to flat geometries. 
For these configurations all the momenta are fixed as functions of the length variables - the 
relations can be obtained by linearizing the formula for the momenta p e valid for flat geometries 
(|8.ip . Furthermore, we have 4-dimensional gauge orbits in the constraint hypersurface. Note 
that these gauge transformations also leave the subspace of flat configurations invariant. Given 
a point p in the subspace of flat configurations there are (N — 4) directions transversal to the 
gauge orbits but tangential to this subspace, i.e. leading to flat configurations which are not in 
the gauge orbit of the point p. There are another (N — 4) directions transversal to the subspace 
of flat configurations but inside the constraint hypersurface. These directions lead to geometries 
with (linearized) curvature. 

These constructions can be enlarged to apply to all constraints at all the vertices of S 
|27j . For the counting of gauge invariant variables note that not all constraints are linearly 
independent. If one considers, for example, the boundary of a 4-simplex as a 3d hypersurface 
£ one counts four constraints at each of the five vertices, hence 20 constraints. Only ten of 
these are linearly independent, however, since the other ten generate global translations and 4d 
rotations of the simplex in the embedding 4d flat space. As there are only ten edge variables in 
S the physical phase space is zero-dimensional. 

10.1 Example: five— valent symmetry— reduced vertex 

Here we consider again the five-valent symmetry-reduced vertex from section 18.11 Using an 
auxilary construction of the 3d star(u) via the gluing of two 4-simplices a(v 1234) and a(v 1235), 
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we already observed in section (|8.1|) that the exterior angle V ; A(i23) & t the triangle A(123) is 
invariant under displacements of the vertex v, that is 



Thus, 



is an observable of the linearized theory. It is the unique observable linear in the configuration 
variables y a ,y h (modulo rescaling). The dihedral angle V ; A(i23) would be the only part of the 
deficit angles in the bulk which depends on the lengths adjacent to the vertex v. Hence, the 
linearized observables are related to the (linearized) deficit angles and ultimately to the identities 
(|4.3|) . expressing the invariance of the deficit angles under the vector fields Yj evaluated on a 
flat background. This invariance, for instance for the deficit angles hinging at the tent pole, can 
be confirmed numerically in this example. 
Explicitly, Va(123) is given by 



2^/6a2 



Y>A(123) = arcsec ^ 3q2 _ 3b 2 + 1 J > < 1 ,K 

so that 



V* = 3V3a(a 2 + b 2 -l) ^ A 

" V-3a* - 36 4 + 6a 2 6 2 + 5a* + 2(^-3 ^ (3a 2 - 1) " " j ' ' 

The Poisson bracket of y^ with the constraint C in (|8. 14[) can be explicitly computed and is 
vanishing. 

Next, we will construct the matrix T and in this way also obtain a momentum observable. 
One choice for T is to define (assuming the generic case det(iV) 7^ 0) 



N gauge a := Y a N aa + Y h ' N ba N gauge 6 := Y a N ab + Y b N bb (10.10) 



and 



/ Y a Y b \ 

Tr e = I * N Y , (10.11) 

\ JV gauge b Jv gauge a J 

where the indices take values T = {gauge, obs} and e = {a, b}. It clearly satisfies the conditions 
in (I10.2p . The inverse is then proportional to 

(T-\ T =( N ^ a ~ Y a b ) . (10.12) 

\ - iv gauge b 1 J 

The observable y obs = J2 e y e {T~ 1 ) e ° bs is proportional to y'^ as defined by (|10.7|) . The 
constraint (ignoring boundary variables) can now be expressed as 

C = Ngauge ~\~ Ngauge gauge V ^ (10.13) 

where N gauge gauge = J2 e , e > Y e N ee ,Y e '. 

The explicit expressions are quite lengthy, but can be computed in a straightforward way. 
For instance, for the specific configuration a = 1,6 = 1 we obtain 

y ^ = vf {~\ ya+2yh ) ' (10 - 14) 
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while the momentum observable, defined via (|1U,1|) . reads 



TT obs = -0.548142 ir a - 0.85287 vr b . (10.15) 
Finally, the constraint ()10.13p is 

C = 7T gauge - 0.294509 y 9auge = . (10.16) 

11 The dynamics of gravitons as generated by the tent moves 

The constraints (|6.5p are linear in the perturbation variables. Besides restricting the allowed 
dynamical configuration, the constraints also generate the changes in the edge lengths and 
momenta induced by infinitesimal displacements of the vertices. However, they leave observables 
- corresponding to the graviton degrees of freedom - invariant, i.e. do not generate any dynamics 
for these. In the continuum the dynamics of the linearized theory with respect to the background 
time (i.e. the time defined by the flat solution and a 3 + 1 decomposition of this solution) is 
generated by a quadratic, global Hamiltonian. Here the dynamics is described by tent moves. 
Also for the discretized case this dynamics is defined with respect to the time as defined by 
the background solution: We evolve the configuration in discrete steps and the proper distance 
between the vertices in different time steps is essentially determined by the background solution. 
The tent move dynamics for the perturbation variables is given by the equations ()6.3[) 

tt* = -M: e ,yi +1 -N: e ,yi-N: b y b , (11.1) 



where the matrices appearing in (jll.ip have been defined in ([6~ 

The momenta at the next time step 7Tg +1 ,7r£ l+1 are defined by (|3.31 13.61) . the linearization of 
which gives 

7r n+1 i, e ' M n 4- N' n i, e ' 4- N' n n b 

n e — Vn M e'e + iV ee'Vn+l + iV ebV 



IT. 



J+i = ^ + A« + N%y* +1 + Nfry" . (11.2) 



Here we additionally introduced the matrices 

d 2 s n d 2 s n ( d 2 s n \ - 1 d 2 s n 



N ,n 



and 



dii +l dii +l oi e n+1 dt n \dt n dt n J dt n dii +l 

d 2 s n d 2 s n ( d 2 s n y 1 d 2 s n 
eb di e n+1 di b di e n+1 dt n \dt n dt n J dt n di b [ • ' 



Nl 



d 2 S n 


d 2 S n ( 


d 2 S n y 


1 d 2 S n 


dl b dl c n 


dl b dt n V 


dt n dt n ) 


dt n dl n 


d 2 S n 


d 2 S n 


( d 2 S n ' 


y 1 d 2 S n 


dl*dl* n+1 


dl b dt n 


\dt n dt n 


) dt n dl e n+l 


d 2 S n 


d 2 S n t 


' o 2 s n y 


- 1 d 2 S n 


dl b dl h ' 


dl b dt n \ 


dt n dt n ) 


dt n dl b ' ' 



W 8l b dl b ' 8l b dt n \dtndtn) 8t n dl b ' ' { ' ' 

The equations (jll.ip have to be used to determine the variables y n+1 as a function of the 
momenta 7r™ and the variables y n , y n - The variables Vn+i a t time step (n + 1) can, however, not 
be uniquely determined as M" , is not invertible and as shown in section [5] has four right null 
vectors n+1 Yf. (Here we introduced an additional index (n + 1) as these are the null vectors 
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defined by the background geometry at time step (n + 1). The left null vectors of Af", are 
n Yf which we so far denoted by just Yf.) This non-uniqueness reflects the gauge freedom in 
the evolution. We will use the splitting of the variables into gauge variant and gauge invariant 
variables introduced in (jlO.ip 

< = "(T- 1 )^ + n (T _1 )g7Tj 

e _ n+l T e a , n+ljie I 
Vn+1 — 1 aUn+l t L lVn+\ 

7T e = (T ) e 1T a + (T ) e TTj . (11.5) 

Equation (jll.ip relating the momenta ir™ and the length variables becomes 

TUrp— X\a„n , n/rp—l\I n -\jrn fn+lrpe'n , n+lrpe',I 

(T ) e 7r a + [T ) e n I = -M ee , T a y n+1 + Tj y. 



ee' ^ a iJn+1 ± I f/n+1 

^ fetf + "Tf'yi) - K b y b ■ (n.e) 



Multiplying this equation with "T^ and remembering that Tf = "Yf and n+1 T| = n + 1 y/ we 
obtain 

vr", = — (^T^,M" e , n+1 T^ y% + i — (^T^N^ e , n T^j y% — ("T^N^) y b , (11.7) 

where we used the conditions (|10.2p on the transformation matrix T. Now the null vectors of 
M™ e , are projected out and we can invert equation (jll.7p for the invariant combinations of the 
length variables y^+x- The y^+i are left undetermined and, consequently, can be freely chosen. 
These variables correspond therefore to lapse and shift. Similarly to equation (111.7p . one can 
show 

<t l = y^^M^^^ + ^^N^^y^i + i^^N'^j/'. (11.8) 
To this end one has to confirm that 

n+v T*,N' n ee , n+x Yf' = (11.9) 
which follows from the condition (|10.2p on the transformation matrix n+1 T 

n+x T e a ,N^ x n+l Yf' = (11.10) 

and the fact that N^ 1 n+1 Yf' = -N' n ee , n+1 l} e ', derived in equation Q . 

With (jll.7p and (jll.8p we reduced the dynamics onto the gauge invariant variables. 

In the rest of this section we will show that the constraints at the next time step are auto- 
matically satisfied. To show the preservation of constraints at the vertex v itself we contract the 
equation (|11.2|) for the momenta 7T™ +1 

vr e " +1 = yiMX + N' n ee ,yi +l + N'ly b (11.11) 



with the vector fields n+1 Yf. Using again N™J~ ln+1 Yf' = -N'l e , n+1 Yf from equation flO} 
and iV™! 1 ™+ 1 y/' = -N' n ee , n + l Yf' from equation §M) we obtain 

n+1 Yf^ +1 = - n+1 YfK e t l y e n + i- n+1 Y I e Kb +1 y b > ( 1L12 ) 

that is the constraints (|6.5|) at time step (n + l). 

We also have to show that the constraints at the neighbouring vertices v' are satisfied after a 
tent move at v has been performed. This is slightly more involved, however, quite straightforward 
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Figure 5: Schematic illustration of the alternating tent moves at neighbouring vertices v and v' , starting 
from the Cauchy hypersurface £ n — l- 

if one starts from the covariant picture. To this end, consider a specific neighbouring vertex v' 
and the situation schematically represented in figure [5j That is, at a time (n — 1) we firstly 
perform a tent move at v' n _ x , then a tent move at v n = u n _i and, finally, at v' n = v' n+l . The 
action for these three tent moves is S = S n -i + S n + S n+ i associated to a piece of triangulation T 
with boundary and one bulk vertex v' n = v' n+l . With S we denote the action with the tent pole 
edges integrated out. To simplify notation we use E as an index for all edges in this triangulation 
(including the boundary) except for the edges e n = e(v' n v n ) and e n+ i = e{v' n v n+ i). With E' we 
denote all edges adjacent to v' n = v' n+1 except for the edges e n = e(v' n v n ) and e n+ \ = e{v' n v n+ i). 
The indices A, A' run over all edges of the triangulation T '■ 

We have one inner vertex in the triangulation and hence four null vectors v Yj for the Hessian 
whose components are v 'Yf', v 'Yf n , «'y/™ +1 . Note that 

" y i JWJI+ Y i j^Ja+ Y i j^TJa (11.13) 

holds for all edges A in the triangulation, also for edges in the boundary. To see this either 
apply arguments similar to those in section [5] or, alternatively, realize that the v Yj would also 
be null vectors for the Hessian associated to a bigger triangulation containing the triangulation 
T ■ But the first derivative in (| 1 1 . 13[) is with respect to edges adjacent to v' n , hence even for the 
bigger triangulation there is only the action S associated to T involved. 

To obtain the constraints at v' n , we firstly contract the equations (|11.13p with the variables 
y A . Secondly, we use the definition for the momenta at v' n 

E , d 2 5 n _i A - en a 2 s w _i A 

and thirdly the equation 
v'ye n+1 _ v'yE' 9 2 S ( d 2 S \ _v'ye- ^ 2 S ( d 2 S 




QlE'Qle n+1 y d le n+1 Qle n+1 I I Qle n Qle n+1 I ^e„ +1 ^e„ +1 

between the components of v 'Yj that follows from the equations (|11.13p with A set to e n+ \. We 
obtain 

= - , Y I E 'Tr E ,+ v 'Yf-7r- en + 



v 'Yj E 'Nl^ +1 y* n+1 + "'if «AL n '" +1 , (11.16) 
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where N AA ^ is given by 

n,n+l = , d 2 (S n + S n+1 ) d 2 S ( d 2 ~S V 1 d 2 ~S 

' A V 1 l A l A ' dl A dl^ I a^+ia/^+i J dl^dl A ' ' 

and coincides with the Hessian of the 'effective action' obtained by integrating out the length 
associated to the edge e n +i from S n + S n +l- This explains also the vanishing of Nj£^~ 1 in the 
last line of (TTJ5| , Hence, these equations do not depend on the variable y en+1 and one can 
check that (| 1 1 . 16 j) indeed reduce to the constraints at v' n at time n. To this end, note that 
according to fjll . 13 j) 

v'vE' AT n,n+l v ' v e n AT n,n+l _ y' v E' ^"Tt-l v ' v e„ 9^Sn-l /ii io\ 

which shows that the constraint equations (111. 16j) . indeed, involve only variables from time step 
n. 

Next, we want to obtain the constraints at the vertex v' n+1 at time step (n + 1) starting from 
the constraints at time step n. We use the defining equations for the dynamics of the tentmove 

CEUIIE2J 

e' e' , d 2 S n A g d 2 S n A d 2 S n A . 

<+i = * n + j^jjy , 7T "+ 1 = -^xv , * n = -pr^y ( n - 19 ) 

and, furthermore, the relation (j 1 1 . 1 5 j) between the components of the null vectors v lj in 'back- 
ward direction'. In the resulting expression all terms involving " Yj n or y e " cancel each other 
and we obtain 

v 'Yt E ' N n+1 v E + v> Y^ n+1 N n+1 v E + v 'y e ' N n+1 7/ 5n+1 4- v 'v^ n+1 N n if n+l C 11 201 

where 

»i* = ws* ■ <"- 21 > 

The equations ()11.2Q|) . in fact, constitute the constraints at v' n+1 . Therefore, the tent move 
dynamics at the vertex v preserves the constraints also at the neighbouring vertices v' . 

The (commutation) algebra of tent moves at neighbouring vertices will be considered in 
further work \27\ 129] , as well as the generalization of a dynamics defined by tent moves to a 
dynamics defined by Pachner moves [30J. 

11.1 Example: symmetry— reduced five— valent vertex 

Using the splitting into gauge invariant and gauge variant variables, the dynamics completely 
decouples. At each time step we have a constraint 

(1 = TT n -4- N n ?/ 9 aM 9 e (-[ 1 OO^i 

^ gauge ' gauge gauge an V 11 '^^/ 

fixing the gauge momentum as a multiple of the gauge variable y n au9e , which, on the other hand, 
can be freely chosen. The gauge invariant variables at different times are coupled through 

''obs Ivl obs obs im+1 ly obs obs an 

^obs 1 = Vn S ^obs obs + ^' obs obs Vn+l (11.23) 
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where M: 



obs obs 



Eee> T obsM2e' T oL ^d SO on. 



Specifically evaluating the last equation for the first two time steps on the flat background 
with a = b = 1 and to = 1/10 yields 



12 Higher order dynamics and pseudo constraints 

In the previous sections we discussed the theory defined by an expansion of the action to 
quadratic order on the flat background. Note that this background solution also displays gauge 
symmetries, namely translations of the vertices. For the linearized theory we obtained con- 
straints - arising as equations of motion which only depend on the dynamical data associated to 
one time step. Although the constraints could have been dependent on the background gauge 
parameters at the next or previous time steps (which do not belong to the dynamical variables), 
they actually did not. In fact, if the constraints were dependent on the background gauge pa- 
rameters at the other time steps, it would have been impossible to obtain consistency of the 
constraint evolution (assuming local evolution laws). 

However, we expect gauge symmetries to be broken for the full non-linear theory [10]. Hence, 
starting with some higher order of the expansion of the Regge action, the gauge freedom should 
become fixed by the equations of motion. Here the lowest order to become fixed is actually 
not the first order variables but the background gauge itself. Interestingly, this fixing of the 
background gauge minimizes the dependence of the Hamilton principal function (sometimes 
called the Hamilton-Jacobi function), i.e. the action evaluated on the solution as a function of 
the boundary data, on the background gauge. 

In the remainder of this section we will discuss this mechanism, which will turn the con- 
straints into pseudo constraints, i.e. equations of motion which depend on lapse and shift |12|, 1 14j . 

Firstly, we will discuss the covariant formulation. The canonical description can be obtained 
afterwards as a rewriting of the equations of motion. In order to simplify the formulae, we will 
ignore variations y b of the edges in the boundary of the tent moves. We consider two consecutive 
tent moves from time step n = to time step n = 2 and consider a boundary value problem 
with data given for times n = 0, 2 and free variables at time n = 1. Moreover, we will assume 
that the lengths of the tent pole edges to,t\ have been integrated out, that is we will work with 
the effective actions So,S\. 

We expand the length variables as l e n = H e n + e ^'y„ + e 2 ^ 2 'y^ + 0(e 3 ) and proceed similarly 
for the momenta. Furthermore, we will use the split of the variables into gauge variant and 
gauge invariant ones defined by the linearized theory, i.e. we will use the transformation matrix 
Tp satisfying the conditions in f) 1 . 2 1) . 

The equations of motion (contracted with x Tp), 



< 6s = -44. 1584 < s - 36.3356 2/0°' 
TT l obs = 44.1584?^ + 31. 0145 yf s 



obs 



(11.24) 









expanded to second order are given by 
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In the following we will consider the equations from the variation of the gauge variables, i.e. 
equations with index T = I. For these the first line in (|12.2p vanishes as it contains the Hessian 
of the action contracted with the null vector Yf = Tf. Denoting the second order terms by &j 
we will show the following 

Claim: &i coincides with the derivative of Hamilton's principal function truncated to sec- 
ond order in the direction of the null vectors 1 Yj. 

Proof: The proof will proceed in two main steps. First we will show that - if we use the 
linearized equations of motion - all terms with gauge variables in 6/ vanish. Consequently, 
there are no variables left to solve for and we have to use equation (|12.2p as a consistency equa- 
tion for the background gauge at n = 1. In a second step, we will show that &i coincides with 
the derivative of the second order Hamilton's principal function with respect to the background 
gauge. 

To begin with, we introduce the notation 

It A' A" \ — 1 d 3 (,go + Si) x , , // « 2 (11, A' (IV A" / 19 o\ 

A / jA » ol l OL n' OL n" 

where A', A" can stand for the gauge indices I" or for the obsevable indices a', a" . 
We start by showing that 

(a) all terms (I\I' n i /\" n „) with n' / 1 and all terms (IiA',/",,) with n" ^ 1 vanish. 

Consider, for instance, (IiIqA",,). The third derivative appearing in this term can be rewritten 



MS 



dlfdl%'dl% 11 1 dl\\ 1 dl e ' dl% 



where we used that % e ,' can be expressed as functions of the length variables /q only (and lengths 
in the boundary of the tent). For n' = 2 the expression in (|12.4p vanishes. For n' = 0, 1 we can 
understand the right hand expression as the double derivative of the momentum at time n' = 
or at time n' = 1, respectively. Since both derivatives are in flat directions in configuration 
space, we can use the expression for 

d So 

Pe"\flat = ("I)™ "^77 \flat (I 2 - 5 ) 

which is valid on flat configurations. On this subspace of the configuration space p™,, is a 
function of the variables 1^ only, hence either the derivative with respect to Zg or the derivative 
with respect to l\ will force the expression in (|12.4p to vanish. 

In conclusion, all second order terms in the equations of motion associated to the gauge 
index I which contain gauge variables at times n = or n = 2 vanish. Note that these terms 
would also vanish if we considered the second order momenta. 

Next, we show that 

(b) all terms (ii/{A"„) and all terms (iiA^,/{') vanish if one uses the first order equations 
of motions for the ^yf . 
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We use a similar rewriting as in (a), that is 

fl»<4 + ft) ^ = a / _ ^8 % A v^sa . (12 . 6) 



n',n"=0,l,2 

1. order eom 



dlfdlfdl% 11 1 dlf \ 1 dlf dl% J V 9l t J dlf'dl^', 



The first term on the right hand side vanishes for the same reason as before: the term inside the 
bracket is zero on flat configurations and the entire expression is a derivative in flat direction of 
this term. Concerning the second term, note that 

E d 2 (S + gi) (i) e " sr^ 9 2 (S + Si) n'^'du" sr^ <9 2 (S + Si) n , a „ 
aie'aie" Vri ~ Z> J " Vn ' + 2^ f)l e ' f)l e " a " Vn ' ^ ' 

n'=0,l,2 m i°V n'=0,l,2 ° l l ° l n' n'=0,l,2 ° l l ° l n' 

are the first order equations of motion associated to the edge e' with the first term on the right 
hand side vanishing automatically. Hence, (ii/{A"„) vanishes if the first order equations of 
motion for the ^yf are satisfied. 

With (a) and (b) the remaining terms in ©/ are then given by 

6/ - (ha' l'{) + (/i/HO + ihlWi) + {hdj'O + (hli<4) + (ha' 2 l'{) 

("=0,1,2 

J2 &<><») ■ ( 12 - 8 ) 
n',n"=0,l,2 

The first line of (|12.8p can be rearranged according to (|12.6p and (|12.7|) to yield terms propor- 
tional to the first order equations of motion for the variables ^'yf ■ If these first order equations 
are satisfied we therefore only remain with terms without any dependence on the first order 
gauge variables and without any second order (gauge and gauge invariant) variables. 

(c) We will consider the second order of Hamilton's principal function - that is the action 
evaluated on the solution - and its derivative with respect to the background gauge parameter. 

The action S = So + Si expanded to second order reads 

S = (0) ^ + £ E ft" 7 ? ( {1 k r + e Vyt) + M % (< V + s <V 

n=0,2 n 1 

+ ^E^"W\ r : ( \ r : + 0(e 3 ) . (12.9) 

n',n" 0i n' al n" 

The zeroth order term does not depend on any (background) variables at time step n = 1, since 
only extrinsic curvature angles appear in it (which can be expressed using length variables from 
only the boundary or time steps n = 0, 2). The same holds for the second term in the first line. 
The last term in the first line vanishes because of the zeroth order equations of motion and also 
its derivative l Yj-^ vanishes even though we would like to solve for the first and second order 
variables. We remain with the second order terms. Using, as in (|12.3p . the notation 

(A ' a" \ ._ i v 92 S n 'r e ' n "r e " p 2 ( x )-7/ a ' Q-XA" (19 ml 

K^n'^n") ■— 2 Z_> file' ft] e" ' J A' 1 A" e Vn' Vn" i (LZ.W) 
A' A" n' OL n" 
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we can analyze the terms according to their type. Firstly, notice that all terms with n' = 
and n" = 2 or vice versa vanish. Secondly, all terms of the type (I{ A"„) and (A' n ,I") vanish 
as l Tf = 1 Yf is a null vector of the Hessian of the action. These term still vanish if we apply 
another derivative ^f-J^ corresponding to infinitesimally changing the vertex at n = 1 in the 

embedding flat space time, where the Hessian contracted with l Yf identically vanishes. The 
same holds for terms of the type {I' n iOi'^„) which vanish either because of (|1U.2|) or ([5. lip . 
We are left with the following second order terms 

= + £ . (12.11) 

n',n"=0,l,2 

The first two terms disappear under the action of a derivative x Yf^ as is shown in (a). For 
the other terms we obtain 

E «'<") = E Ui«>"») + e • (12-12) 

1 n',n"=0,l,2 n',n"=0,l,2 

The additional terms summarized as <B arise through the derivative ^Yf^s acting on the solu- 
tions for ,^?/f and on the components 1 T^,, 1 T^„ . (We have to replace the variables 
^ ( l )y& by the solutions to obtain Hamilton's principal function. Also note that the deriva- 
tives with respect to the length l\ are not acting on the components n T^, for n' ^ 1. The reason 
is that the expression n YfN" e , only involves background variables from time step n' . Similarly, 
the conditions (|10.2p on the matrix n Tp only involve background variables from time step n', 
hence one can also choose n Tp to be of this type.) These terms £ are proportional to the first 
order equations of motion, however, and therefore vanish 

& - V 2 V d2s n 'r e ' l r e " W?/"' ( l Y e d ( l \, a "\ 4- 

d 2 s 



12 ^ 8 S l , n n ,, A 8 (l) a >\ (1) a " 

2 £ 2^ WW' 1 ot' 1 a" [ Y I fl]e Vl j Vn" + 
n"=0,\,2 m ^ m n" V Ml J 

ST 32 S n> T e> ( lye A ^T e "\ ^\, a ' i l \, a " 4- 
,^ 1)2 d%dlf ^ { Yl dlf T «") ^ ^ + 



±e 2 

2 fc 



x 2 d S ( x 8 j / » n „ „ (1) / (1) // 

0. (12.13) 



1. order eom 

As a consequence, we finally obtain 



1. order eom 

This finishes the proof. 



61 = l Y[-^ {2) ~S ■ (12.14) 



To summarize, for the first non-linear order of the equations of motion (|12.2p the following 
situation arises: the equations for V = a have to be used to determine the second order gauge 
invariant observables ^y a , as these only appear there. For the remaining equations of motions 
T = I, which contain only first order gauge invariant variables (if the first order equations of 
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motion are satisfied), we do not have any variables left to solve for and we seem to have an 
inconsistent theory. 

However, the remaining terms in the equations of motion for T = I will generically depend 
on the background gauge parameters, which in a sense are zeroth order variables. Indeed, 
these equations of motion now have a precise interpretation, namely as equations which fix the 
background gauge such that the second order of Hamilton's principal function (which can also 
be called the effective action) depends minimally on this gauge. 

This also entails that one can obtain a consistent expansion to higher order only for certain 
choices of the gauge parameter in the background solution. For other choices one cannot expand 
the fluctuation variables y e in a power series in e: the solutions corresponding to the gauge de- 
grees of freedom have a lowest order term proportional to e^ 1 , which can be interpreted as a 
change of the background gauge. On the other hand, at the lowest non-linear order we find that 
the first and second order gauge variables y remain undetermined. For the next order, i.e. an 
expansion of the action to fourth order, we expect that the equation of motion determine the 
first order gauge variables. 

These considerations can be tested with the parametrized harmonic oscillator (and unhar- 
monic generalizations). The action for one time step is given by 

5'n = ^%-^ t4- - \u{q n + q n+ i) 2 (t n+1 - t n ) . (12.15) 

^ \tn+l ~ t n) o 

We consider the variation of the variables at time step n = 1 with fixed data at time steps 
n = 1, 2 and expand the action using q n = e + e 2 ^q n and t n = (°\ n + £^'t n + e 2 ( 2 h n to 
third order around the configurations qo,Qi,Q2 = and arbitrary t^- These configurations are 
solutions to the equations of motion with tk being the background gauge parameters. One finds 
that the second order equation of motion corresponding to t\ is satisfied automatically only for 
t\ = \(t$ + £2)- As a result, the higher order terms determine the time discretization. 

Moreover, if one just defines q^ = + ey^ and tj, = + £Z k an d takes the expansion of 
the action to third order as a definition of the dynamics one finds that the solution of z is not 
analytic in e. The lowest order rather scales with e _1 , that is effectively changes 

In the canonical framework we can define the momenta at time step n = 1 via the action 
Sq and the action S\. The contraction of these momenta with the null vectors Yf resulted 
in constraints. From the previous discussion we can conclude that the second order momenta 
(defined via Si) contracted with ^Yf are of the form 
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Note that only gauge variables ^y 1 , ^y 1 from time step n = 1 appear. The only variables 
from time step n = 2 are the first order gauge invariant ^Vf ■ Using the first order equations of 
motion (jll.6p . however, these can be expressed as linear combinations of variables ^yf, 
at time step n = 1. 
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Hence, if we consider only the fluctuation variables with k > 1 as true variables 

we can also obtain at second order relations which only involve the variables at one time step. 
From this point of view one can still speak of constraints. However, these constraints are 
not automatically preserved by time evolution anymore. The reason is that the corresponding 
covariant equations (|12.8p . which are exactly the condition for the preservation of the constraints, 
are not automatically satisfied. 

If we also consider the (gauge) parameters of the background solution as zeroth order vari- 
ables, however, the second order terms of the constraints (|12. 16[) will in generic cases depend 
on these variables from time step n = 2. In this sense the second order constraints are pseudo 
constraints. Not all the terms in (112. 16D will depend on the background gauge parameters at 
time step n = 2 - one can show that all terms with first or higher order gauge variables only 
depend on the background variables at time n = 1. These are exactly the terms which cancel 
automatically in the covariant equations of motion f|12.8j) . 

Using a similar rewriting as in equation fjl2.6f) and the first order equations of motion, the 
second order part of the constraints can be written as 

, ( lye 9 lye'\ f^-Ua" (1) V (1)1 ( lye 9 wA ^ S l lye' (1)1' (1) a " 



l d 3 Si x , e , , e > (1) p (1) p, 
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lV 19 ^PP? 1 a ' 11 Vn> Vn» , (12-17) 

n ,n"=l,2 

where appearing in the last line can be substituted by an expression involving only variables 
at n = 1 

lrpe d 2 <Sl 2rre' (l\,a' 1 l T e_^l W (l)„,a' /i, lfi \ 

T «dqdif % ' 2/2 = ^ ~ Ta difdi? Ta ' 2/1 • (12 - 18) 

Note that the constraints might remain exact constraints, i.e. relations between variables 
(including zero order variables) from only one time step, to higher or even all orders. The latter 
is the case for tent moves at four-valent vertices, which lead to flat dynamics. The full non-linear 
constraints for this situation are given by (|7.1ip . 

The phenomenon that not all solutions of the linearized theory can be completed to solutions 
of the full theory is similar to the occurrence of linearization instabilities in continuum general 
relativity for space times with compact spatial slices [28J. In our case this phenomenon occurs 
because the solutions of the full theory are unique (if there are no flat vertices), whereas the 
linearized solutions admit firstly freedom for the choice of the background gauge and secondly 
freedom in the choice of the first-order gauge parameters. The consistency conditions eliminate 
this gauge freedom order by order. There is an important difference to the linearization instabil- 
ities appearing in [28J: whereas there the additional conditions are on the first order (physical) 
modes, the consistency conditions here fix the zeroth-order background gauge variables (and 
presumably the higher order equations fix the higher order gauge variables). 

12.1 Example: Symmetry reduced five valent vertex 

Here we will consider the consistency equation arising from the second order equations of motion 
for the symmetry reduced five valent vertex, described in section (|8.ip . The variables we have 
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to deal with are the lengths a n ,b n for n = 0,1,2. We fix data at n = 0, 2 and consider the 
equations of motion with respect to ai, b\. 

To begin with, we have to (numerically) find solutions for the lengths of the tent pole to, t\. 
This can be done to second order in an expansion around the flat configuration determined by 
the initial values ^clq = 1, (°'bo = 1 and 

{ % = ^+r, ( % = i~r (12.19) 

where r is the background gauge parameter, determining the position of the vertex at n = 1 
in the background flat space time. Note that ^02,^62, the background data at n = 2 are 
independent of r. 

This way, we can obtain the effective action S expanded to third order. From this action 
we can obtain the first and second order equations of motion for ^yf, ^~'y\, ^yf, ^y\- Solving 
the first order equation corresponding to the derivative with respect to a\ for ^yf and using 
the solution in the other first order equation of motion, one will find that it is automatically 
satisfied. This is the signature of the exact gauge freedom for the linearized theory. 

Using the first order solution we can solve the second order equation corresponding to the 
derivative with respect to a\ for ( 2 'yf. Again, we use this solution in the other second order 
equation corresponding to the derivative with respect to b\. For instance, for r = we find that 
we have to solve the equation (ignoring terms of order 1CP 8 arising due to numerical errors) 

= -0.5464921 ( (1) ?/o) 2 - 0.9715415 ( (1) ?/o) 2 + 1-4573123 ( (1) yo)( (1 Vo) 
+1.0936310 ( (1) y a )( (1) y 2 a ) - 1-3269178 (W^X^) 
-1.4581747 ( (1) ^)( (1) y 2 a ) + 1.7692237 (C^X* 1 ^) 

-0.5471391 ( (1) y 2 ) 2 - 0.8054603 {^ ] y b 2 ) 2 + 1.3277030 ( {1) y 2 )( (1) y 2 ) • (12.20) 

However, as expected from the previous discussion, all variables at time step n = 1 have dropped 
out. We thus have to find a value for the background gauge parameter r such that the remaining 
second order equation is satisfied. A priori one would expect that the value of r has to depend on 
the boundary data ^Vn> ^Un with n = 0, 2. But (as for the parametrized harmonic oscillator) it 
turns out, that this equation can be solved independently from these first order boundary data. 
All coefficients in (|12.20p which are non-vanishing for r = vanish simultaneously (within 
numerical accuracy) for r = —0.008303982. We conjecture that there is a general mechanism 
which ensures that the second order terms 6/ can be made to vanish independently of the value 
of the first order boundary data. 

13 Conclusion and discussion 

In this work we introduced a canonical formalism for discretized gravity which exactly repro- 
duces the dynamics as defined by the discrete action. In this way also the exact and approximate 
symmetries of the action are reproduced as first class constraints and pseudo constraints, re- 
spectively. For linearized Regge calculus, which exhibits exact gauge symmetries, we obtained 
Abelian constraints valid for arbitrary triangulated Cauchy surfaces. The momenta and con- 
straints are local functions (as opposed to the suggestion in |12j) and have an immediate geo- 
metric interpretation as generators of vertex displacements. The constraints can be shown to 
commute and to be preserved by the tent move dynamics as defined by the quadratic order of 
the action. Although this can be proved by staying entirely in the canonical framework, it is 
much easier to remember that the constraints follow from the symmetries of the action. For 
instance, the constraints are Abelian because second derivatives (of the action) commute. 

Discretized constraints are often derived by performing first a canonical analysis of the 
continuum action and then discretizing the resulting constraints. This often leads to a change 



33 



of the constraint algebra from first to second class [13, [T2], i.e. the constraints close only modulo 
terms proportional to some power in the lattice spacing. This work suggests to not only qualify 
discretized constraints under this criterium but to also consider the properties of the constraints 
near solutions relevant for the continuum limit - here the flat solutions. 

We argued that the constraints to the first non-linear order will aquire some dependence on 
the background gauge parameters and turn therefore into pseudo constraints [2] . These pseudo 
constraints are not automatically preserved under tent move evolution anymore. Rather some of 
the equations of motion turn into consistency conditions, selecting a specific background gauge 
and ensuring that the constraints are now preserved under time evolution. 

These consistency conditions can be rewritten as the derivatives with respect to the back- 
ground gauge parameters of the quadratic part of Hamilton's principal function. Note that this 
quadratic part can be computed in the linearized theory. Despite the fact that the linearized 
theory displays exact constraints and gauge symmetries (in the fluctuation variables), the as- 
sociated Hamilton's principal function is not independent of the background gauge, or in other 
words the discretization. These findings entail that a consistent perturbative expansion is only 
supported for specific background choices. For other choices solutions are not analytic in the 
expansion parameter. 

Although the background gauge is fixed by the lowest non-linear order dynamics the first 
and second order gauge variables remain undetermined. We expect that these are fixed one after 
the other in the higher order dynamics, so that in the n-th order dynamics the nth and (n— l)th 
order of the gauge variables remain undetermined. 

We want to underline that although we found pseudo constraints starting with second order 
for Regge calculus, there might be alternative discretization schemes which lead to first class 
discrete constraints for the full theory. In fact, as was shown in j2] for Regge calculus with a 
cosmological constant, where the symmetries are broken, it is possible to find an alternative 
discretization of the action with exact symmetries also for curved solutions. The methods used 
in [2] are the so-called perfect actions, which exactly reproduce the dynamics (and therefore 
symmetries) of the continuum. These perfect actions can be obtained by a coarse graining proce- 
dure (from the continuum). Obtaining the perfect action for 4d gravity will be very complicated, 
but might be achievable in a perturbative approach [29]. The considerations in this paper give 
necessary prerequisites towards this end. 

If one is only interested in a particular class of solutions, for instance almost homogeneous 
solutions for cosmological applications, one might try to derive improved actions by adapting 
the discretization to the chosen background. In the case of cosmology, it would be interesting 
to have an action which displays gauge symmetries for homogeneous solutions and in this way 
ensure the correct number of propagating degrees of freedom on cosmological backgrounds. See 
also [31] for related discussions in loop quantum cosmology and |32] for a proposal for using tent 
moves in a (quantum) cosmological setting. 

We hope that these results can be useful for connecting covariant and canonical approaches 
to quantum gravity, in particular spin foam models and canonical loop quantum gravity. This 
has been achieved in 3d [33] (where the symmetries are exact for the full non-linear theory) but 
is an outstanding problem in 4d [34 j . Spin foam models can be understood as partition functions 
for a discretized theory [35]. In canonical loop quantum gravity the central ingredient defining 
the dynamics is given by the quantum constraints |36j . 

This work suggests a derivation of a canonical quantum theory from the covariant one, that is 
a canonical theory directly derived from the amplitudes associated to spin foam models. As spin 
foam models are based on a discretization of the Plebanski action one expects the symmetries to 
be broken. Hence, one can doubt whether on the discrete level, i.e. before taking any continuum 
limit or sum over triangulations, (exact) quantum constraints exist for the full theory that would 
reproduce the dynamics defined by the spin foam model. Nevertheless, it should be possible to 



34 



define a canonical dynamics with discrete time that reproduces the spin foam amplitudes. To 
this end it would be useful to perform a canonical analysis of the discretized Plebanski action, a 
task that was only recently achieved for the continuum action in |37j . In particular, in |17j one 
can find the relation between the canonical phase spaces of discretized Plebanski theory and of 
Regge calculus. Also in [17] it was pointed out that there exists a class of triangulations which 
admits only flat solutions and on which the symmetries are therefore exact. For this class of 
triangulations it might therefore be possible to derive quantum constraints from the quantum 
amplitudes. To accomplish this, the symmetries of the amplitudes have to be better understood. 
Interesting results in this direction and a connection to the quantum constraints can be found 
in [381 [39]. 
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